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Abstract: We compute the supersymmetric (SUSY) contributions to the observ- 
ables in B° s —>■ K + K~ and B° s — > K°K° decays. The hadronic parameters in the 
standard-model (SM) amplitudes are obtained from the B® — > K K decay using a 
recent approach that combines flavor SU(3) symmetry and a controlled input from 
QCD factorization. The latest experimental data for BR(B Q S — > K + K~) is in agree- 
ment with the SM prediction. We study how the branching ratios and the direct and 
mixing-induced CP asymmetries of both B° s — > KK decay modes are affected with 
the inclusion of SUSY, after imposing constraints from BR(B — > Xg^), B — > ttK 
and AM S over the parameter space. While the branching ratios remain unaffected 
by SUSY, we identify the CP asymmetries of the B® — > KK decays as the most 
promising observables to look for large deviations from the SM. 
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1. Introduction 

Recently, the precision of -B-physics measurements has increased dramatically due to 
the experimental results of CDF, Babar and Belle. This implies that new strategies 
are necessary for controlling hadronic uncertainties. In addition, it is important 
to identify those observables which are useful for signalling the presence of physics 
beyond the standard model (SM). Once these are found, the next step is to explore 
the impact of well-motivated models. 

In a recent publication |] we computed the supersymmetric (SUSY) contribu- 
tions to B° s — > K + K~ decays. In the present paper, we present an update of this 
analysis, with four important improvements. First, we extend the calculation to in- 
clude Bg — > K°K° decays. Second, our discussion of the predictions of the SM is 
based on a new method which uses the B% — > K°K° decay 0. It combines QCD 
factorization with flavour symmetries, and represents a substantial improvement in 
the control of hadronic uncertainties. Third, the limits on the SUSY parameter space 
include the latest constraints from B° s -B° s mixing |J, as well as data from B —>■ Y s 7 
and B — > ttK decays. The fourth point is related to squark mixing. For each fermion, 
each of the two components (left-handed, right-handed) has a scalar SUSY partner, 
the squark. One can have mixing of the left-handed or right-handed squarks of the 
three generations. In Ref. 0, it was assumed that one of the two mixings (LL or 
RR) was zero. In this paper, we allow simultaneous nonzero values of both LL and 
RR mixing. 

We pay particular attention here to the CP-violating asymmetries of the decays 
B° s -> K+K- and B° s -> K°K°. These are defined as: 

\A\ 2 -\A\ 2 lm(e- l ^A*A) , . 

-fiair ._.»>, -fornix " r, . — ,9 j V / 

1-41 +\A \A+\A 
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where A is the amplitude of the decay in question, A is the amplitude of the CP- 
conjugate process, and (fi s is the phase of B®-B® mixing. 

In the SM, B° s -> K°K° is dominated by a single penguin decay amplitude. 
There are other penguin contributions, but they are suppressed and enter only at 
the level of < 5%. (They are all included in Sec. 2.) The direct CP-violating 
asymmetry, A<nr, involves the interference of the dominant penguin amplitude and 
the suppressed contributions, and is therefore very small in the SM. As for the mixing- 
induced CP asymmetry, A m i X , since the B° s — > K°K° decay is dominated by a single 
amplitude, A mix essentially measures <p s . This phase is very small in the SM [Q, so 
that this CP asymmetry is expected to be correspondingly small. Since both CP 
asymmetries of B® — > K°K° are expected to be so small in the SM, this makes them 
interesting observables for detecting the presence of new physics (NP). The situation 
is somewhat different for the decay B° s — > K + K~ since it receives a tree contribution 
which cannot be neglected with respect to the dominant penguin contribution. In 
the SM, the interference of the penguin and tree amplitudes in B° s — > K + K~ gives 
rise to larger CP asymmetries than in B° s — > K°K°. Both decays involve a b — > s 
transition and therefore have branching ratios of O(10 -5 ). 

All of these SM predictions can change in the presence of SUSY. Naively, one 
would guess that all SUSY contributions to B® — ► K + K~ and B° s — > K°K° are 
suppressed by M^/Mg USY , where M SUSY ~ 1 TeV, and are therefore small. However, 
some of the SUSY contributions involve squark-gluino loops, which are proportional 
to the strong coupling constant a s . Compared to the SM, the relative size of these 
contributions is therefore (a s /a)(M^/Mg USY )- Since this is 0(1), such contributions 
can compete with those of the SM, leading to significant modifications of the SM 
predictions for B® — > K + K~, and especially for B® — > K°K°. In this paper, we 
consider only these SUSY contributions, as they are the dominant effects. 

Experimentally, the branching ratio of B® — > K + K~ has been measured at 
CDF[|: 

BR(B° S -> K + K-) cxp = (24.4 ± 1.4 ± 4.6) x 1(T 6 . (1.2) 

The CP asymmetries for B° s — > K + K~ are likely to be measured soon at CDF |J. 
However, the measurements for B® — > K°K° will probably take more time at CDF 
or may have to wait until LHCb. 

In Sec. 2, we discuss the SM expectations for the various observables in B° s — > 
K + K~ and B° s — > K°K°. In order to do this, we must consider certain B% decays. 
Here we follow Ref. || and use B® — > K°K°. Sec. 3 contains the general analysis 
of B° s -> K + K~ and B° s -> K°K° decays with the addition of NP. We discuss the 
amplitudes for these decays in the presence of NP, as well as strategies for measuring 
the NP parameters. We turn specifically to SUSY in Sec. 4 and calculate its effect 
on the amplitudes of B° s -> K+R- and B° s -> K°K°. We note that SUSY can 
substantially modify the SM predictions for the CP asymmetries while keeping the 
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branching ratios basically unaffected. We conclude in Sec. 5. 



2. Standard-Model Analysis of B° s K + K~ and B° s K°ii: 

Consider first B° s — > K + K~ , which at the quark level is 6 — ^ suu. This decay receives 
a contribution from a penguin diagram pen' (the prime indicates a b — > s transition). 
The penguin diagram receives contributions from each of the internal quarks u, c and 
t. However, using the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
one can write 

pen' = V: b V us {P'u - PI) + V; b V cs (Pc - PI) • (2.1) 

The decay B° s — > K + K~ also receives several other diagrammatic contributions [7|]. 
The most important of these is the tree amplitude, which is proportional to V* b V us : 
tree' = V* b V us T'. The amplitude can therefore be written 1 

A(B° S -> K+K-) ~ V: b V us [T' + {P' u - P t ')] + V: b V cs {P' c - P/) 

= V: b V us T s± + V; b V cs P s± , (2.2) 

where P s± = {P' c - P[) and T s± = [T' + (i* - P/)]. Note that |V„* 6 K S | ^ 5%| V^V^I, 
and this CKM suppression compensates the relative size of the amplitudes \P S± /T S± \ ~ 
0.1 (see Ref. 0). Thus, the first term is smaller than the second, but must be in- 
cluded in the analysis. 

The amplitude for B° s —>■ K°K° (quark level: b — > sdd) can be treated similarly. 
In this case, there is no tree diagram, but we keep the notation T s0 for the penguin 
contribution proportional to V* b V us , leading to 

A(B° S -> K°K°) ~ V: b V us T s0 + V; b V cs P s0 . (2.3) 

Here, both T s0 and P s0 are of the same size, but since nothing compensates for the 
strong CKM suppression, the first term is strongly suppressed, leading to a very 
small direct CP asymmetry. 

In the isospin limit, P s± = P s0 . However, T s± ^ T s0 due to the presence of the 
tree diagram in B° s — > K + K~ . Since the terms proportional to V* b V us are small, the 
amplitudes for B° s — > K + K~ and B° s —>■ K°K° are approximately equal in the SM. 
However, this need not be the case for NP. 

In order to determine the amplitudes for B° s — > K + K~ and B° s — ► K°K°, we need 
to know P s± , P s0 , T s± and T s0 . These can be obtained by considering P° decays. 
One can use P° — > 7r + 7r~ decays in order to determine the hadronic parameters 
M, |, 0, [TO, O. However, it was noted in Ref. § that B° d -> if°i? provides 



1 We use the following notation (see Ref. Q): quantities carrying the superscripts dO, sO and s± 
correspond to -> S° -> il /? and £?° -> K+R-, respectively. 
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smaller errors on the predictions, and that the inclusion of some (quite generic) sign 
predictions from B d — > 7r + 7r~ decays leads to a greater restriction on the SM ranges. 

The argument leading to the determination of the B d — > K°K° parameters is as 
follows. The amplitude for this decay can be written 

A(B° d - K°K°) ~ V^KdT* + V^P* ■ (2.4) 

There are thus three unknown quantities to be determined: the magnitudes of P d0 
and T d0 , and their relative strong phase. Three pieces of information are therefore 
needed. One piece of information comes from the measurement of the B d —>■ K°K° 
branching ratio: BR(B° d -> K°K°) = (0.96 ± 0.25) x 10~ 6 [|TJ 2 . 



A second piece of information comes from QCD factorization (QCDf) fl4] , [T5 ] . 
In QCDf the various hadronic quantities can be calculated using a systematic ex- 
pansion in l/m6. However, a potential problem arises when one encounters endpoint 
infrared (IR) divergences in higher-order power-suppressed terms. Their evaluation 
thus requires an arbitrary IR cutoff, and they may be enhanced numerically (for a 
given IR cutoff). The key observation || is that the difference A^ = T d0 — P d0 is 
free of these IR divergences and can be calculated fairly accurately within QCDf: 

A d = (1.09 ± 0.43) x 10~ 7 + i(-3.02 ± 0.97) x 10~ 7 GeV . (2.5) 

Note that the values of the real and imaginary pieces of A^ can be affected by a global 
phase transformation, so that only the modulus |A^| is physical. This provides the 
second constraint on the hadronic parameters of B d — > K°K°. 

Finally, the authors of Ref. find that only values —0.2 < A d dir < 0.2 are 
consistent with the measured value of BR(B d — > K°K°) and the theoretical value of 
A^. This is the third piece of information. 

Using the values of the branching ratio and |A<f|, as well as the allowed range for 
A^ r , one can obtain the moduli and relative strong phase of the hadronic parameters 
in B° d -> K°K°: 

\T d0 \ = (1.1 ± 0.8) x 10" 6 GeV , \p d0 /T d0 \ = 1.2 ± 0.2 , 

arg(P d0 /T d0 ) = (-1.6 ± 6.5)° . (2.6) 

In fact, there is a twofold discrete ambiguity in determining these quantities, 
but the authors of Ref. argue that only one solution is physical. The argument 
uses the two methods involving the decays B d — > K°K° and B d — > 7r + 7r~. It was 
shown in Ref. @ that the second solution (the unphysical one) requires a large U- 
spin violation (in the phase) and, moreover, it predicts A s d f r < 0. This is clearly in 



2 While completing this work, this measurement has been updated by BABAR in Ref. |Tj| to 
BR(B° d -> K°K°) = (1.08 ± 0.30) x 10~ 6 . However, in this work we prefer to stick to the quoted 
value more near to the HFAG value from ICHEP06 
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contradiction with the prediction for the sign of AS; using B® — > as can be 

seen in Table 1 of Ref. |lj] . This resolves the two- fold ambiguity. 

In addition, it was pointed out that there exists a strong anticorrelation between 
the signs of A s mix and A^ r: as can be seen in Table 1 of Ref. fl2| . On the other hand, 
using the method of B® — > ix + n~ one can see that the predicted value for A s ^ ix is 
always negative Jl|. Putting together both methods, one finds an important restric- 
tion on A^ r that prefers positive values. This affects and improves all predictions, 
shown in Table 1 of Ref. [[| (only the lower half of this table should be taken). 

Recently, a measurement was reported by the BABAR collaboration [Tj 

A dir (B d -> K° S K° S ) = -0.40 ± 0.41 (2.7) 

This very preliminary measurement, although still quite uncertain, seems to point 
towards negative values for A^, but it is compatible with a small positive asymmetry. 
Note, however, that the measured central value for A^ ix is much larger than 1. Thus, 
one should really take these numbers only as a proof that they can be measured. 
Here we present the results for the range -0.2 < Af ir < +0.2 (this includes the 
non-preferred negative region). 

From these, one can now compute the parameters of the B® — > K°K° amplitude, 
taking into account SU(3) breaking [0. Factorizable SU(3)-breaking corrections are 
introduced by means of P s0 = fP d0 and T s0 = fT d0 , where the factor / is defined 
as ||15 

Ml o F B °- R \0) 
f = ' B o^ RO = 0-94 ± 0.20 . (2.8) 
M 2 B0 F C ?<- K (0) 



d 



and the input values are taken from Refs. j|, [Uj. This parameter can be calculated on 
the lattice. Other sources of SU(3) breaking which are suppressed by 1/m^ originate 
from hard-spectator scattering (differences in the distribution amplitudes of B^ and 
B® ) and weak annihilation (diagrams in which the gluon emission comes from the 
spectator quark). All effects are computed within QCDf (we assume that QCDf gives 
at least the right order of magnitude), which gives the following bounds @: 



P s0 /(fP d0 ) - 1| < 3% 



o 



|T s 7(/r dU ) - 1| < 3% . (2.9) 

Note that large final-state- interaction SU(3)-breaking effects are possible. However, 
these are common to both B® s —>■ K°K° decays and, consequently, they cancel in 
relating the two modes. Thus the parameters of the amplitude for B® — > K°K° can be 
established (with errors), and the branching ratio and CP asymmetries determined. 
These numbers are given below. 

The decay B® — > K + K~ is somewhat more complicated, as a combination of U- 
spin and isospin is required to connect B^ — > K°K° to B° s — > K + K~ . The relations 
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between those hadronic parameters are given in Ref. ||, including SU(3)-breaking 
corrections evaluated within QCDf. The bounds are 

\P s± /(fP d0 ) - 1| < 2% , 

\T s± /(A s KK ai) ~ 1 -T d0 /(^«i)l < 4% , (2.10) 

where A S KK and 6t\ are additional hadronic parameters that can be estimated in 
QCDf. However, it is also noted that since the T s± term is CKM suppressed, any 
uncertainties in the determination of A S KK and a\ affect the branching ratio and CP 
asymmetries of B° s — > K + K~ only marginally. 

Putting all this together, the SM predictions for the branching ratios and CP 
asymmetries in B° s — > K + K~ and B° s — > K°K° can be obtained. If one conservatively 
takes all values of A^ r between —0.2 and 0.2, the prediction for the branching ratio 
is | : 

BR(B° S -> K°K°) = (18 ± 7 ± 4 ± 2) x 10~ 6 , 
BR(B° S -> K + K-) = (20 ± 8 ± 4 ± 2) x 10~ 6 . (2.11) 

But if only values of A d ^ ir > (up to 0.2) are taken, the prediction becomes: 

BR(B° S -> K°K°) = (18 ± 7 ± 4 ± 2) x 10" 6 , 
BR(B° S K + K~) = (17 ± 6 ± 3 ± 2) x HT 6 . (2.12) 

The first error reflects the uncertainty in the QCDf estimates of and «i, as well as 
in BR(B® — > K°K°) (this is the largest uncertainty). The second error corresponds 
to the uncertainty in / (SU(3) breaking). The third error introduces a rough estimate 
of non-enhanced l/m^-suppressed contributions. 

The predicted ranges for the CP asymmetries of B° s — > K + K~~ and B° s —>■ K°K° 
within the SM are illustrated in Fig.|l|. The conservative predictions are [[|: 

-0.011 < Af ir < 0.011 , 

-0.015 < A s £ ix < 0.005 , 

-0.22 < Af r < 0.49 , 

-0.55 < A*, < 0.40 . (2.13) 

Or, considering only positive values of A^.: 

-0.011 < Af ir < 0.003 , 
-0.015 < AZ X < 0.005 , 
-0.22 < A s d f r < 0.49 , 

-0.55 < A s ± x < 0.02 . (2.14) 

As expected, the CP asymmetries for B® — > K°K° are predicted to be very small in 
the SM, with the prediction that \Af ir \ should be equal or less than 1% and l-A^I 
should be less than 2%. 
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Figure 1: SM predictions for the CP asymmetries in B® — ► K°K° (up) and B® — » K + K 
(down) as a function of -Adir(-B^ — * K°K°). As explained in the text, the preferred range 
is the non-shadowed half of the plots [A dir (B^ -» i-T°K°) > 0]. 



3. New Physics 

At present, there are many measurements of B decays and several discrepancies with 
the SM have appeared. For example, the CP asymmetry in b — > sqq modes (q = 
u, d, s) is found to differ from that in b — > ccs decays by 2.6a (they are expected to be 



approximately equal in the SM ) fig , [17| . In addition, some B — ► nK measurements 



disagree with SM expectations fl8| , although the so-called B — > txK puzzle |TTJ has 



been reduced pOfl . One also sees a discrepancy with the SM in triple-product 
asymmetries in B — > <pK* pl| , |22| , and in the polarization measurements of B — > 
0if* |23j, g§ and 5 — >■ pif* [F2£| 0. Although these discrepancies are not yet 
statistically significant, there is a unifying similarity: they all point to new physics 
in b — > s transitions. We will therefore follow this indication and assume that the 
NP appears in b — > s decays but does not affect b — > d decays. That is, B® — > K + K~ 
and B° s -> K°K° can be influenced by the NP, but B G d -> K°K° |2| is not. 



There are many NP operators which can contribute to a given B decay. How- 
ever, in Ref. [29|, it was observed that the matrix elements of NP operators carrying 



non-negligible strong phases are necessarily suppressed with respect to the SM con- 
tribution. (Note that each NP contribution can in principle have a different strong 
phase.) As a result, the relevant NP matrix elements can be combined, and a given 
B decay thus receives a single NP contribution. 
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The amplitude for the decay B® — ► K + K can therefore be written 

A(B° S -> K+iT-) = ^± + .4V*" . (3.1) 

In the above, A s ^ contains both weak and strong phases, but we have separated out 
the weak phase (<& u ) of the NP contribution. (Note: the NP strong phase is zero, so 
that A u is real.) The name of the NP parameters reflects the fact that this decay is 
b — > suu at the quark level. 

Similarly, the amplitude for B° s — > K°K° (quark level: b — ► sdc?) is 

-4(fi° -> K°^°) = 4f M + . (3.2) 

If the NP conserves isospin, we have A u = A d and $ u = $d, but in general this need 
not be the case. 

One can make experimental measurements of B® — > K + K~, obtaining the 
branching ratio and the direct and mixing-induced CP asymmetries. Assuming that 
A s s ^j is known from B® — > K°K° decays (as in Sec. 2), these three measurements 
allow one to extract A u and $ u , as well as the relative strong phase between the SM 
and NP amplitudes. One can proceed similarly for B® —>■ K°K°. In this way one 
can measure all the NP parameters 0. 

In the following section, we compute the SUSY contributions to A u , A d and 
and the resultant branching ratios and CP asymmetries for B® — > K + K~ and 
B° S ^K°K°. 



4. SUSY Predictions for B° K+K" and -> K°ii: 

As mentioned above, the relevant SUSY contributions to 6 — > sgg transitions come 
from squark-gluino box and penguin diagrams. We follow the procedure outlined in 
Ref. [|T[, which is based on the work of Grossman, Neubert and Kagan PP |. We refer 
the reader to these references for further details. 

As shown in Ref. |PH| , the most natural solution to the B — > ttK puzzle is the 



introduction of isospin-breaking NP amplitudes. The isospin-breaking effect is more 
naturally realized in the present scenario, which allows large up-down squark mass 
splittings, than in the more popular mass insertion (MI) approximation. Indeed, since 
we are working in a scenario with near-maximal mixing between bottom and strange 
squarks, the squark MI approximation is not adequate. Note that the dangerous 
s — > d and b — > d flavour-changing neutral currents are not generated in this scenario 
due to the assumption of vanishing (1, 2) and (1, 3) components in the scalar down- 



type squark mass matrix |30] 



In this scenario, the SUSY contributions to the Wilson coefficients depend on 
the following parameters: the masses of the squarks and gluino, two mixing angles 
9l,r and two weak phases Sl,r- These angles parametrize the rotation matrices that 
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diagonalize the left- and right-handed squark mass matrices. The expressions for 
the NP amplitudes A u e^ u and A d e q ' d in terms of these parameters are obtained in 
complete analogy with Ref. ||J]. To be specific, we write the expressions for these 
amplitudes in terms of the Wilson coefficients (denoted by cf and Cg? ): 

Gp 



.4 



V2 

where x — 1-18 and 

A = (^ |(fod)v + A|-B°)(ir°|(Js)v + A|0) = %{m\ - m 2 K )f K F B ^ K ~ 1.42 GeV 3 (4.2) 
For the numerical values we take |, |S| = 0.31 and / K = 0.1598 GeV. The 



explicit expressions for the Wilson coefficients cf ,u (which includes both L and R 
mixing contributions) can be found in the appendix of Ref. [ffl, where some small 
typos in Ref. |3(J were corrected. 

In order to obtain the NP amplitudes, we must evaluate the Wilson coefficients 
within SUSY. We consider the following values and ranges for the SUSY parameters: 

• rriu L = m- dhR = m~ HR = m~ g = 250 GeV 

• 250 GeV < m^m- SRL < 1000 GeV 

• -7T < 5 LiR < 7T 

• -tt/4 < 6 L)R < it/4 

The SM inputs are the same as those used in Ref. 0. The Wilson coefficients are 
sensitive to the s — b mass splitting. They vanish for m s - = m~ b and grow when 
the splitting is large. We therefore expect these contributions to be most important 
for large values of m s - (keeping m~ b fixed). In the same way, NP effects in b — > dqq 
transitions depend on the difference mj — m~ b . By setting m d - = m~ b we ensure that 
b — ► d decays get no such contributions, which is consistent with the discussion in 
Sec. 3. A difference between A u e® u and A d e' s?d is only possible in the presence of 
a nonzero u — d mass splitting. Without it there are no contributions to isospin- 
violating operators. However, this mass splitting must be very small in the left- 
handed sector due to SU{2) L invariance. We therefore set = mj L , but allow for 
a significant mass splitting in the right-handed sector. 

There are also constraints on the SUSY parameter space from other processes 
that have been already measured. The constraints from the decays B — > ttK and 
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B — > X s j are described in Ref. Jl|. In particular, we take BR(B — > X^) 



(3.55 ± 0.26) x 10 -4 [16|, with an increased error to cover the various theoretical un- 
certainties 3 . Most importantly, one must take into account the recent measurement 
of AM S , which was not included in the previous analysis. The latest measurement 



together with the SM fit |33|, |35| , gives 



(0.81 ± 0.19) ps" 1 (4.3) 



AM| M 



exp 



The constraints from all these measurements have been included in our analysis. 
Other traditionally-important constraints like B® — > are very sensitive to other 

SUSY parameters, mostly tan/3 and uia- However, for small values of tan/? and 
values for uia above 200 GeV, they have no effect on our allowed region to SUSY 
(Fig. 0). 

Taking into account the various constraints, the contributions from LL and RR 
mixing have been analyzed. AM S is the strongest constraint, and it is the relevant 
one when considering only LL or RR mixing separately. In particular, it has a large 
impact on the phases $ u and In the case of B° s — > K + K~ , LL mixing gives the 
largest contribution to the amplitude, more than twice that of RR. However, in the 
case of B® — > K°K° both contributions are similar in size. 

When both LL and RR mixings are allowed simultaneously, the constraints on 
the SUSY parameter space are changed. In this case new operators for B° s -B° s mixing 
are generated, so that the effect is not simply the combination of the two separate 
contributions (for instance, see Ref. [0). We find that (i) now BR(B — > Xg'j) is 
also important, not only AM S , and (ii) the global effect of the constraints is weaker. 
The upshot is that there is a certain enhancement of the NP amplitudes when both 
LL and RR mixings are combined. In particular, the weak phases $ u and <3>d are not 
so strongly constrained as when either LL or RR mixing is taken to vanish. 

In Fig. |2| we show the allowed ranges for A u e % ^ u and A d e l ® d in the scenario 
with simultaneous LL and RR mixings. The dark regions correspond to the values 
that these amplitudes take when varying the parameter space over the initial ranges. 
The light regions show how these values are reduced by the existing experimental 
constraints mentioned above. There are two important remarks. First, we see that 
the above constraints do indeed greatly reduce the allowed SUSY parameter space. 
Second, even so, the effect on A u e l ® u and A d e % ^ d can be significant. 

At this stage we can identify what are the effects of the various constraints in 
reducing the SUSY parameter space. The bound from B — > X s 7 affects only the 



3 Interestingly, the latest NNLO calculations in the SM show that BR(B -> X s j) (SM) is a little 
lower than the experimental average pS2f . 

4 Note that here we take the largest fit result for the SM prediction ^] because it falls within 
the second prediction of Ref. Moreover, the value for AAf| M taken here differs from that used 
for the bound in Ref. [fl| that was based on Ref. [B6|- 
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Figure 2: SUSY contribution to the NP amplitudes ,4V*" (left) and .A<V* d (right) in 
the scenario with simultaneous LL and RR mixings. The dark regions correspond to the 
variation of the SUSY parameters over the considered parameter space. The light regions 
satisfy the experimental bounds, including the recent measurement of AM S . 

left-handed sector. In particular, for large rag, the regions with \9l\ > 10°, \5l\ < 60° 
and \9l\ > 10°, \Sl — vt | < 60° are excluded. The bound from AM S is much stronger: 
when rris > 400 GeV, any values of \9l,r\ > 5° are excluded, as well as those regions 
in which 5l + Sr ~ — 37r/2, —n/2, it/2, 37r/2. 5 After these bounds are imposed on the 
parameter space, the constraints from B — > ttK have very little effect on the regions 
in Fig. 0. 

Note that the allowed region for A u is much larger than that for A d , by approx- 
imately a factor of 3. In the isospin limit, these should be equal, so this factor of 3 is 
a measure of isospin breaking in this NP scenario. In particular, for = 250 GeV 
(zero UR-d,R mass splitting), the values of A u reduce to those for A d . 

We now examine the effect of these contributions on the observables. By adding 
the SUSY contributions to the SM amplitudes as in Eqs. ( |3. 1| ) and ( |3.2| ), it is pos- 
sible to compute the branching ratio and the CP asymmetries in the presence of 
SUSY. Fig. |3] shows the allowed values for the B° s — ► K + K~ observables, for three 
different values of Af^ r , compared with the predictions of the SM and with the re- 
cent experimental value for the B° s — ► K + K~ branching ratio reported by the CDF 
collaboration || [Eq. (|1.2j) 1. The agreement between the CDF measurement and the 
prediction of the SM in Ref. erases any discrepancy between experiment and the 
SM. This branching ratio will now be an important future constraint. The branch- 
ing ratio within SUSY should not deviate much from the SM prediction so as not 

5 For further details of the AM S constraint on this parameter space, see Ref. p4| . 
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K+K-) - A diT (B° s -» K+K-) (up) and A mix (B° s -» - Aiir(£° if+if") (down) 

in the presence of SUSY, for a) A% = -0.1, (b) = and (c) Ag r = 0.1. The dashed 
rectangles correspond to the SM predictions. The horizontal band shows the experimental 
value for BR(B° S -► K + K~) at la. 



to generate any disagreement with data. Indeed, Fig. |3] shows that the impact of 
SUSY on the branching ratio of B° s — > K + K~ is practically negligible. Interestingly, 
for positive values of A^ r (preferred region), the SM predicts a smaller value for 
BR(B° S — > K + K~), but it is now compatible with the new data. Still, it is in this 
case that SUSY shows a larger deviation in the correct direction. 

A completely different picture arises for the CP asymmetries. The results for 
the direct CP asymmetry reveal that SUSY can have an impact. This is not surpris- 
ing: SUSY introduces a term in the total amplitude which is of the same order of 
magnitude as that of the SM and carries a weak phase that is not constrained. The 
mixing-induced CP asymmetry gets affected in a more dramatic way. The interpre- 
tation is that the SUSY contribution to the mixing angle cj) s can be large (in fact it 
can take all values between —n and 7r), while in the SM it is tiny: 0f M ~ —2°. Any 
experimental measurement falling inside the dark area in the plots, but outside the 
dashed rectangle, would not only signal NP but clearly could be accommodated by 
sup er symmet ry. 

Fig. [| shows the results for B® —>■ K°K°. Although the branching ratio is little 
changed in the presence of SUSY, the enhancement of the CP asymmetries due to 
the inclusion of the SUSY contributions is in this case even more important. The 
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K°K°) - A dir (B° s -> K°K°) (up) and A mix (B° s -» K°if°) - 4^(5° -> ^°^°) (down) in 
the presence of SUSY, for (a) Af ir = -0.1, (b) = and (c) = 0.1. The dashed 
rectangles correspond to the SM predictions. These are quite small in the three lower plots, 
so they are indicated by a circle. 



reason is that, within the SM, the CP asymmetries are much smaller in B° s — ► K°K° 
than they are for B® — > K + K~ , because of the absence of the tree diagram. Thus 
the impact of SUSY is much greater. This is evident by looking at the three lower 
plots in Fig. |4], where the tiny rectangles corresponding to the SM predictions can 
hardly be observed. We have drawn a circle around them to indicate their position. 

These are a good illustration of the general scenario discussed in the introduction. 
While the branching ratios in this case are relatively insensitive to supersymmetry, 
the direct and mixing-induced CP asymmetries of these decays are greatly affected. 
Thus, these CP asymmetries are the observables to focus on in order to observe 
NP, particularly SUSY, while the branching ratio of B® — > K + K~ can become an 
important constraint on models beyond the SM other than SUSY. 



5. Conclusions 

In this paper we have considered the branching ratios and CP asymmetries for the 
decays B° s — > K + K~ and B° s — > K°K° in a supersymmetric (SUSY) model, focusing 
on the dominant gluino-squark contributions | 30fl . This analysis is an extension and 



update of Ref. ffl, which now includes the SUSY analysis of both B® decay modes 
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and allows for both LL and RR mixing. In addition, the determination of the SM 
contributions from the decay B® — > 7r + 7r~ used in Ref. [|J has been replaced by the 
recently-proposed combination of B® — > K°K° and QCD factorization of Ref. J|. 
This allow us to obtain more precise predictions. 

We have included the constraints coming from BR(B — > X s j), B — ► 7rK and 
AM S , and we find the following results. 

• The new-physics (NP) amplitudes are A u e % ® u (B® — > K + K~) and ^4 d e** d 
(5° — > K°K°). We find that both can get significant contributions from SUSY. 
In the isospin limit, these quantities are equal. However, our calculations show 
that, for the region of parameters considered, in SUSY there can be a difference 
of up to a factor of 3 between the NP amplitudes. This indicates the possible 
level of isospin breaking in this type of theory. In particular, in the SUSY 
model considered here, large isospin violation is possible when there is large 
mass splitting in ur-cIr. 

• The branching ratio BR(B° S — ► K + K~) is very little affected by SUSY. At 
most, the SM prediction can be increased by 15% for A^ r = 0.1. In fact, 
SUSY can somewhat improve the already good agreement between the SM 
prediction and the new precise CDF measurement ||. The impact of SUSY 
on BR(Bg — > K°K°) is even smaller, reflecting the reduced allowed region for 
A d e % ^ d as compared to A u e l ® u . 

• The situation is very different for the CP asymmetries; the size of the effect 
depends strongly on the decay and the type of asymmetry. For B® — ► K + K~ , 
the direct CP asymmetry within SUSY is in the range —0.1 < A c u r (Bg — > 
K+K-) SUSY < 0.7 for -0.1 < Af ir < 0.1. Depending on the value of Af r , it 
may be possible to disentangle the SUSY contribution from that of the SM. 
This is due to the competition between the tree and the NP amplitudes for 
each value of A^ r . As for A m i X (B° s — > K + K~), its value can vary all the way 
from —1 to +1, signaling a large impact from SUSY. 

• Turning to B° s — > K°K°, the CP asymmetries are particularly promising. This 
decay is dominated by the penguin amplitude in the SM, and so the direct CP 
asymmetry is strongly suppressed: it is predicted to be at most of the order of 
1%. However, in the presence of SUSY, the direct CP asymmetry can be 10 
times larger. The mixing- induced CP asymmetry is also predicted to be very 
small in the SM. However, A mix (B° s -> K°K°) SUSY covers the entire range, 
and so this asymmetry can be large in the presence of SUSY. 
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